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ANALYSIS OF LANDING-GEAR BEHAVIOR ! 


By BENJaMIN Miuwrrzky and Francis E. Соок 


SUMMARY 


This report presents a theoretical study of the behavior of the 
conventional type of oleo-pneumatic landing gear during the 
process of landing impact. The basic analysis is presented in 
a general form and treats the motions of the landing gear prior 
to and subsequent to the beginning of shock-strut deflection. In 
the analysis of the first phase of the impact the landing gear 18 
treated as a single-degree-of-freedom system in order to deter- 
mine the conditions of motion at the instant of initial shock-strut 
deflection, after which instant the landing gear 18 considered as 
a system with two degrees of freedom. The equations for the 
two-degree-of-freedom system consider such factors ав the 
hydraulic (velocity square) resistance of the orifice, the forces 
due to air compression and internal friction in the shock strut, 
the nonlinear force-deflection characteristics of the tire, the wing 
lift, the inclination of the landing gear, and the effects of wheel 
spin-up drag loads. 

The applicability of the analysis to actual landing gears has 
been investigated for the particular case of a vertical landing 
gear in the absence of drag loads by comparing calculated 
results with experimental drop-test data for impacts with and 
without tire bottoming. The calculated behavior of the landing 
gear was found to be in good agreement with the drop-test data. 

Studies have also been made to determine the effects of varia- 
tions in such parameters as the dynamic force-deflection 
characteristics of the tire, the orifice discharge coefficient, and the 
polytropic exponent for the air-compression process, which 
might not be known accurately in practical design problems. 

The study of the effects of variations in the tire characteristics 
indicates that in the case of а normal impact without tire 
bottoming reasonable variations in the force-deflection character- 
istics have only а relatively small effect on the calculated behavior 
of the landing gear. Approximating the rather complicated 
force-deflection characteristics cf the actual tire by simplified 
erponential or linear-segment variations appears to be adequate 
for practical purposes. Tire hysteresis was found to be 
relatively unimportant. In the case of a severe impact involving 
tire bottoming, the use of simplified exponential and linear- 


segment approximations to the actual tire force-deflection 
characteristics, which neglect the effects of tire bottoming, 
although adequate up to the instant of bottoming, fails to indicate 
the pronounced increase in landing-gear load that results from 
bottoming of the tire. The use of exponential and linear- 
segment approximations to the tire characteristics which take 
into account the increased stiffness of the tire which results from 
bottoming, however, yrelds good results. 

The study of the importance of the discharge coefficient of the 
orifice indicates that the magnitude of the discharge coefficient 
has a marked effect on the calculated. behavior of the landing 
gear; a decrease in the discharge coefficient (or the product of the 
discharge coefficient and the net orifice area) resulis 1п an 
approximately proportional increase in the maximum upper- 
mass acceleration. 

The study of the wmportance of the air-compression process 
in the shock strut indicates that the air springing às of only 
minor significance throughout most of the impact and. that 
variations in the effective polytropic exponent n between the 
isothermal value of 1.0 and the near-adabatic value of 1.8 have 
only a secondary effect on the calculated behavior of the landing 
gear. Even the assumption of constant air pressure in the strut 
equal to the initial pressure, that 18, n=0, yields fairly good 
results which may be adequate for many practical purposes. 

In addition to the more exact treatment, an investigation has 
been made to determine the extent to’which the basic equations 
of motion can be simplified and still yield acceptable resulte. 
This study indicates that, for many practical purposes, the 
air-pressure force in the shock strut can be completely neglected, 
the tire force-deflection relationship can be assumed to be linear, 
and the lower or unsprung mass can be taken equal to zero. 
Generalization of the equations of motion for this simplified 
system shows that the behavior of the system 18 completely 
determined by the magnitude of one parameter, namely the 
dimensionless wnitial-velocity parameter. Solutions of these 
generalized equations are presented in terms of dimensionless 
variables for а wide range of landing-gear and impact parameters 
which may be useful for тата йу estwunating ИМЕЙ 
performance in preliminary design. 


1 Supersedes NACA TN 2755, “Analysis of Landing-Gear Behavior’ by Benjamin Milwitzky and Francis E. Cook, 1952. 
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INTRODUCTION 
The shock-absorbing characteristics of airplane landing 


gears are normally developed largely by means of extensive - 


trial-and-error drop testing. The desire to reduce the ex- 
pense and time required by such methods, as well as to pro- 
vide a more rational basis for the prediction of wheel-inertia 
drag loads and dynamic stresses in flexible airframes during 
landing, emphasizes the need for suitable theoretical methods 
for the analysis of landing-gear behavior. Such theoretical 
methods should find application in the design of landing 
gears and complete airplane structures by permitting 

(a) the determination of the behavior of a given landing- 
gear configuration under varying impact conditions (velocity 
at contact, weight, wing lift, etc.) 

(b) the development of a landing-gear configuration to 
obtain a specified behavior under given impact conditions 

(с) a more rational approach to the determination of wheel 
.spin-up and spring-back loads which takes into account the 
shock-absorbing characteristics of the particular landing gear 
under consideration 

(d) improved determination of dynamic loads in flexible 
airplane structures during landing. This problem may be 
treated either by calculating the response of the elastic sys- 
tem to landing-gear forcing functions determined under the 
assumption that the airplane is a rigid body or by the simul- 
taneous solution of the equations of motion for the landing 
gear coupled with the equations representing the additional 
degrees of freedom of the structure. In many cases the former 
approach should be sufficiently accurate, but in some 
instances, particularly when the landing-gear attachment 
points experience large displacements relative to the nodal 
points of the flexible system, the latter approach, which takes 
into account the interaction between the deformations of the 
structure and the landing gear, may be required in order to 
represent the system adequately. 

Since many aspects of the landing-impact problem are so 
intimately connected, with the mechanics of the landing gear, 
the subject of landing-gear bebavior has received analytical 
treatment at various times (see bibliography). Many of 
the earlier investigations, in order to reduce the mathematical 
complexity of the analysis, were limited to consideration of 
highly simplified linear systems which have little relation to 
practical landing gears. Some of the more recent papers 
consider, with different degrees of simplification, more real- 
istic nonlinear systems. The present report represents an 
attempt at a more complete analysis of the mechanics of 
practical landing gears and, in addition, investigates the im- 
portance of the various elements which make up the landing 
gear, as well as the extent to which the system can be reason- 
ably simplified for the purpose of rapid analysis. 

The basic analysis is-presented in a general form and takes 
into account such factors as the hydraulic (velocity square) 
resistance of the orifice, the forces due to air compression and 
internal friction in the shock strut,.the nonlinear force- 
deflection characteristics of the tire, the wing lift, the inclina- 
tion of the landing gear and the effects of wheel spin-up drag 
loads. An evaluation of the applicability of the analysis to 
actual landing gears is presented for the case of a vertical 
landing gear in the absence of drag loads by comparing cal- 
culated results with drop-test data. 
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Since some parameters, such as the dynamic force- 
deflection characteristics of the tire, the orifice discharge 
coefficient, and the polytropic exponent for the air-compression 
process, may not be accurately known in practical design 
problems, a study is made to assess the effects of variations 
in these parameters on the calculated landing-gear behavior. 

Studies are also presented to evaluate the extent to which 
the dynamical system can be simplified without greatly im- 
pairing the validity of the calculated results. In addition to 
the investigations for specific cases, generalized solutions for 
the behavior of a simplified system are presented for a wide 
range of landing-gear and impact parameters which may be 


useful in preliminary design. 
. SYMBOLS 

Ag pneumatic area 

А, hydraulic area 

А, area of opening in orifice plate 

4А, internal cross-sectional area of shock-strut inner 
cylinder 

Ag external cross-sectional area of shock-strut inner 

i cylinder 

A, cross-sectional area of metering pin or rod in 
plane of orifice 

P. TES net orifice area 

Са orifice discharge coefficient 

d overall diameter of tire 

F, . pneumatic force in shock strut 

Е, hydraulic force in shock strut 

Е, friction force in shock strut 

Fs total axial shock-strut force 

F, normal force on upper bearing (attached to inner 
cylinder) 

Р, normal force on lower bearing (attached to outer 
cylinder) | 

Fy, force normal to axis of shock strut, applied at 
axle 

Fy, vertical force, applied at axle 

Ен, horizontal force, applied at axle 

Fg, . resultant force, applied at axle 

Fs, force parallel to axis of shock strut, applied to 
tire at ground 

Fy force normal to axis of shock strut, applied to 


tire at ground 
Fy, vertical force, applied to tire at ground 


Fy, horizontal force, applied to tire at ground 

Ев, resultant force, applied to tire at ground 

g gravitational constant 

КІ lift factor, L/W 

L lift force 

ІІ “ axial distance between upper and lower bearings, 
for fully extended shock strut 

ls axial distance between axle and lower bearing 
(attached to outer cylinder), for fully ox- 
tended shock strut 

a,b,m,r constants corresponding to the various regimes 
of the tire-deflection process 

а” combined constant, ad 

m combined constant, md’ 

n polytropic exponent for air-compression process 


in shock strut 
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R Reynolds number 

Da air pressure in upper chamber of shock strut 

D» hydraulic pressure in lower chamber of shock 
strut 

Q volumetric rate of discharge through orifice 

ra radius of deflected tire 


9 shock-strut axial stroke 

T wheel inertia, torque reaction 

t time after contact 

T time after beginning of ain seat deflection 

v air volume of shock strut 

16 polar moment of inertia for wheel assembly 


about axle 

Vy vertical velocity 

^H horizontal velocity 

W total dropping weight 

W, weight of upper mass above strut 

И». weight of lower mass below strut 

ta horizontal displacement of lower mass from 
position at initial contact 

2 vertical displacement of upper mass from рові- 
tion at initial contact | 

2 vertical displacement of lower mass from ровї- 
tion at initial contact 

Uy dimensionless upper-mass displacement from 
position at initial contact 

Ue dimensionless lower-mass displacement from 
position at initial contact 

т dimensionless shock-strut stroke, u—us 

0 dimensionless time after contact 


angle between shock-strut axis and vertical 


Omas 
f a" ас 
로 


Ns shock-strut effectiveness, ~; 
Ші marmar 


uy 
Су; da 
Nig landing-gear effectiveness, ^, 
6 time interval in numerical integration procedures 
m coefficient of friction. between tire and runway 
Hi coefficient of friction for upper bearing (attached 
to inner cylinder) 
Ua coefficient of friction for lower bearing (attached 
to outer cylinder) 
p mass density of hydraulic fluid 
a angular acceleration of wheel 
Axes: 
г vertical axis, positive downward 
2 horizontal axis, positive rearward 
Subscripts: 
0 at instant of initial contact 
T at instant of initial shock-strut deflection 
su at instant of wheel spin-up 
тах maximum value 
Notation: 
IC 9] absolute value of ( ) 


( )* estimated value of ( ) 

The use of dots over symbols indicates differentiation with 
respect to time ¢ or г. 

Prime marks indicate differentiation with respect to 
dimensionless time 0. 
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MECHANICS OF LANDING GEAR 
А DYNAMICS OF SYSTEM 


In view of the fact that landing-gear performance appears 
to be relatively unaffected by the elastic deformations of 
the airplane structure (see, for example, refs. 1 and 2) par- 
ticularly since in many cases the main gears are located 
fairly close to the nodal points of the fundamental bending 
mode of the wing, that part of the airplane which acts on a 
given gear can generally be considered as a rigid mass. 
As a result, landing-gear drop tests are often conducted in 
а jig where the mass of the airplane is represented by a 
concentrated weight. In particular instances, however, such 
as in thé case of airplanes having large concentrated masses 
disposed in an outboard position in the wings, especially 
airplanes equipped with bicycle landing gear, consideration 
of the interaction between the deformation of the airplane 
structure and the landing gear may be necessary to repre- 
sent the system adequately. 

Since the present report is concerned primarily with- the 
mechanics of the landing gear, it is assumed in the analysis 
that the landing gear is attached to a rigid mass which has 
freedom only in vertical translation. The gear is assumed 
infinitely rigid in bending. The combination of airplane 
and landing gear considered therefore constitutes a system 
having two degrees of freedom (see fig. 1(а)) as defined by 
the vertical displacement of the upper mass and the vertical 
displacement of the lower or unsprung mass, which is also 
the tire deflection. The strut stroke s is determined by 
the difference between the displacements 2, and 2» and, in 
the case of inclined gears, by the angle о between the axis 
of the strut and the vertical. For inclined gears, compression 
of the shock strut produces a horizontal displacement of the 
axle 2ҙ. From consideration of the kinematics of the system 


2;)tan e. 


In the analysis, external lift forces, corresponding to the 
aerodynamic lift, are assumed to act on the system through- 
out the impact. In addition to the vertical forces, arbitrary 
drag loads are considered to act between the tire and the 
ground. 

The system treated in the analysis may therefore be con- 
sidered to represent either a landing-gear drop test in a jig 
where wing lift and drag loads are simulated, or the landing 
impact of a rigid airplane if rotational motions are neglected. 
Rotational freedom of the airplane, where significant, may 
be taken into account approximately by use of an appro- 
priate effective mass in the analysis. 

Figure 1 (b) shows a schematic representation of a typical 
oleo-pneumatic shock strut used in American practice. The 
lower chamber of the strut contains hydraulic fluid and the 
upper chamber contains air under pressure. The outer cyl- 
inder of the strut, which is attached. to the upper mass, 
contains a perforated tube which supports a plate with a 
small orifice, through which the hydraulic fluid is forced to 
flow at high velocity as a result of the telescoping of the 
strut. The hydraulic pressure drop across the orifice thus 
produced resists the closure of the strut, and the turbulence 
created provides a powerful means of absorbing and dis- 
sipating a large part of the impact energy. In some struts 
the orifice area is constant; whereas, in other cases a metering 


е 2 — 2 . 
it can be seen that s— 때 2 and m=s sin o— (21— 
p 
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(a) 
(a) System with two degrees of freedom. 


ж = (4-22) ton Ф 
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(b) 


(b) Schematic representation of shock strut. 


Кісовв 1.—Dynamical system considered in analysis. 


pin or rod is used to control the size of the orifice and govern 
the performance of the strut. 

'The compression of the strut produces an increase in the 
air pressure which also resists the closure of the strut. In 
figure 1 (b) p, represents the oil pressure in the lower chamber 
and ра represents the air pressure in the upper chamber. 

In addition to the hydraulic resistance and air-pressure 
forces, internal bearing friction also contributes forces which 
can appreciably affect the behavior of the strut. 

The forces created within the strut impart an acceleration 
to the upper mass and also produce an acceleration of the 
lower mass and a deflection of the tire. Figure 1 (c) shows 
the balance of forces and reactions for the wheel, the inner 
cvlinder, and the outer cylinder. It is clear that the strut 


and the tire mutually influence the behavior of one another 
and must be considered simultaneously in analyzing thesystem. 
FORCES IN SHOCK STRUT 
From consideration of the pressures acting in the shock 
strut it can be readily seen from figure 1 (b) that the total 
axial force due to hydraulic resistance, air compression, and 
bearing friction can be expressed by 


Ез=р(4,—А,) + pa (da — А) рад, P, 


where 

А internal cross-sectional area of inner cylinder 

Ag external cross-sectional area of inner cylinder 

A, cross-sectional area of metering pin or rod in plane 


of orifice 


~ 
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Forces on outer cylinder 


_ 21-42 
7 COS p 


5 


Forces on inner cylinder 


(с) 


Forces on wheel 


(c) Balance of forces and reactions for landing-gear components. 
Ғтаувв 1.—Concluded. 


-- ~ ----- —- - - -- -- 
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This expression can also be written йз 
Fs= (Pr— pa) (Ài— Ap) + peda +P, 
--(Рь- Ра) Ant Polat Fy 
=F +F +F, : (1) 


where | 

Pn—Pa pressure drop across the orifice 

А, hydraulic area (А; — 4, for the strut shown in fig. 1) 
A. pneumatic area (Аз for the strut shown in fig. 1) 

In this report the terms (p,—,)A, and Paa are referred 
to as hydraulic force F, апа pneumatic force F,, respec- 
tively. For the strut shown-in figure 1, the hydraulic 
and pneumatic areas are related to the strut dimensions as 
previously noted. In the case of struts having different 
internal configurations, the hydraulic and pneumatic areas 
may bear somewhat different relations to the dimensions 
of the strut. In such cases, however, consideration of the 
pressures acting on the various components of the strut 
should permit these areas to be readily defined. 

Hydraulic force.— The hydraulic resistance in the shock 
strut vesults from the pressure difference associated with the 
flow through the orifice. То. a landing gear the orifice area 
is usually small enough in relation to the diameter of the 
strut so that the jet velocities and Reynolds numbers are 
sufficiently large that the flow is fully turbulent. As a 
result the damping force varies as the square of the tele- 
scoping velocity rather than linearly with the velocity. 
Since the hydraulic resistance is the major-component of 
the total shock-strut force, viscous damping cannot be 
reasonably assumed, even though such an assumption 
would greatly simplify the analysis. | 

The hydraulic resistance can be readily derived by making 
use of the well-known equation for the discharge through 


an orifice, namely, 
Q— Cans Z (Pa— Pa) 
where 


Q volumetric rate of discharge 
C, coefficient of discharge 
A, net orifice area 
2, hydraulic pressure in lower chamber 
p, ait pressure in upper chamber 
p шазз density of hydraulic fluid 
From considerations of continuity, the volumetric rate of 
discharge can also be expressed as the product of the tele- 
scoping velocity 8 and the hydraulic area 4, | 


Q= A,’ 


Equating the preceding expressions for the discharge per- 
mits writing the following simple equation for the pressure 
drop across the orifice 
mp pA,’ 8? 
D» Pe 9(0,4,)° 


The hydraulic resistance Е, due to the telescoping of the 
strut is given by the product of the differential pressure 
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Pa—pa and the area A, which is subjected to the hydraulic 
pressure, ав previously noted. Thus 


А». 
f= XC AF 8? (2) 


Equation (2) can be made applicable to both the compres- 
sion and elongation strokes by introducing the factor ET 
to indicate the sign of the hydraulic resistance; thus 


8 Ap . 
FQ—— p 4X5 2 


=$] ЗА 8 (2a) 


The net orifice area A, may be either a constant or, when a 


. metering pin is used, can vary with strut stroke; that is, 


A,=A,—A,=A,(s), where A, is the area of the opening іп 
the orifice plate and A, is the area of the metering pin in the 
plane of the orifice. At the present time there appears to bo 
some tendency to eliminate the metering pin and uso a con- 
stant orifice area, particularly for large airplanes, in which 
case A =Á.. In the general case, the orifice discharge 
coefficient might be expected to vary somewhat during an 
impact because of changes in the size and configuration of 
the net orifice area, changes in the exit conditions on the 
downstream face of the orifice due to variations in the amount 


- of hydraulic fluid above the orifice plate, changes in the entry 


conditions due to variations in the length of the flow chamber 
upstream of the orifice, and because of variations in the 
Reynolds number of the flow, so that, in general, C,=0,(8, R). 
Although the individual effects of these factors on the dis- 
charge coefficients for orifices in shock struts have not been 
evaluated, there is some experimental evidence to indicate 
appreciable variations of the discharge coefficient during 
impact, particularly in the case of struts with metering pins. 
It might be expected that such variations would be con- 
siderably smaller for gears having a constant orifice area. 

In order to evaluate the precision with which the orifico 
discharge coefficient has to be known, a brief study is 
presented in a subsequent section which shows the effect of 
the discharge coefficient on the calculated behavior of a 
Janding gear with а, constant orifice area, under the assump- 
tion that the discharge coefficient is constant during the 
impact. 

The foregoing discussion has been concerned primarily with 
the compression stroke of the shock strut. Most struts 
incorporate some form of pressure-operated rebound check 
valve, sometimes called a snubber valve, which comes into 
action after the maximum stroke has been attained and closes 
off the main orifice as soon as the strut begins to elongate, so 
that the fluid is forced to return to the lower chamber through 
small passages. The action of the snubber valve introduces 
greatly increased hydraulic resistance to dissipate the energy 
stored in the strut in the form of air pressure and to prevent 
excessive rebound. The product C,A, to be used in equation 
(2a) during the elongation stroke is generally uncertain. The 
exact area A, during elongation is usually somewhat difficult 
to define from the geometry of the strut since in many cases 
the number of connecting passages varies with stroke and the 


- 


leakage area around the piston may be of the same order of 
magnitude as the area of the return passages. Furthermore, 


the magnitude of the orifice discharge coefficient, and even. 


possibly the nature of the resistance, are questionable due to 
the foaming state of the returning fluid. Fortunately, the 
primary interest is in the compression process r&ther than 
the elongation process since the maximum load always occurs 
before the maximum strut stroke is reached. 

Pneumatic force.—The air-pressure force in the upper 
chamber is determined by the initial strut inflation pres- 
sure, the area subjected to the air pressure (pneumatic area), 
and the instantaneous compression ratio in accordance 
with the polytropic law for compression of gases, namely 


pave" == Constant, or 
рег. (>) 
where 


Pa іг pressure in upper chamber of shock strut 

Pao air pressure in upper chamber for fully extended strut 
v air volume of shock strut 

7% X air volume for fully extended strut 

Since the instantaneous air volume is:equal to the difference 
between the initial air volume and the product of the stroke 


and pneumatic area Áa, p,— Pa ( 해 2 
ji 


to the air pressure is simply the product of the pressure and 
the pneumatic area: 


Еа-тра 4 C ў (3) 


In the preceding equations, the effective polytropic 
exponent n depends on the rate of compression and the rate 


of heat transfer from the air to the surrounding environment. 


Low rates of compression would be expected to result in 
values of n approaching the isothermal value of 1.0; whereas 
higher values of n, limited by the adiabatic value of 1.4, 
would be expected for higher rates of compression. The 
actual thermodynamic process is complicated by the violent 
mixing of the highly turbulent efflux-of hydraulic fluid and 
the air in the upper chamber during impact. On the one 
hand, the dissipation of energy in the production of turbu- 
lence generates heat; on the other hand, heat is absorbed by 
the aeration and vaporization of the fluid. The effect of this 
mixing phenomenon on the polytropic exponent or on the 
equivalent air volume is not clear. A limited amount of 
experimental data obtained in drop tests (refs. 3 and 4), 
however, indicates that the effective polytropic exponent 
may be in the neighborhood of 1.1 for practical cases. A 
brief study of the importance of the air-compression process 
and the effects which different values of n may have on the 
calculated behavior of the landing gear is presented in a 
subsequent section. 

Internal friction force.—In the literature on machine de- 
sign the wide range of conditions under which frictional 
resistance can occur between sliding surfaces is generally 
classified in three major categories, namely, friction between 
dry surfaces, friction between imperfectly lubricated surfaces, 
and friction between perfectly lubricated surfaces. In the 
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case of dry friction, the resistance depends on the physical 
characteristics of the sliding surfaces, is essentially propor- 
tional to the normal force, and is approximately independent 
of the surface area. The coefficient of friction и, defined as 
the ratio of the frictional resistance to the normal force, is 
generally somewhat greater under conditions of rest (static 
friction) than under conditions of sliding (kinetic friction). 
Although the coefficient of kinetic friction generally de- 
creases slightly with increasing velocity, it js usually con- 
sidered, in first approximation, to be independent of velocity. 
If, on the other hand, the surfaces are completely separated 
by a fluid film of lubricant, perfect lubrication is said to exist. 
Under these conditions the resistance to relative motion 
depends primarily on the magnitude of the relative velocity, 
the physical characteristics of the lubricant, the area, and 
the film thickness, and is essentially independent of the 
normal force and the characteristics of the sliding surfaces. 
Perfect lubrication is rarely found in practice but is most 
likely under conditions of high velocity and relatively small 
normal pressure, where the shape of the sliding surfaces is 
conducive to the generation of fluid pressure by hydro- 
dynamic action. In most practical applications involving 
lubrication, a state of imperfect lubrication exists and the 
resistance phenomenon is intermediate between that of dry 
friction and perfect lubrication. 

In the case of landing-gear shock struts, the conditions 
under which internal friction is of concern usually involve 
relatively high normal pressures and relatively small sliding 
velocities. Moreover, the usual types of hydraulic fluid 
used in shock struts have rather poor lubricating properties, 
and the shape of the bearing surfaces is generally not con- 
ducive to the generation of hydrodynamic pressures. It 
would therefore appear that the lubrication of shock strut 
bearings is, at best, imperfect; in fact, the conditions appear 
to approach closely those for dry friction. In the present. 
analysis, therefore, it is assumed, in first approximation, that 
the internal friction between the bearings and the cylinder 
walls follows laws similar to those for dry friction; that is, 
the friction force is given by the product of the normal force 
and a suitably chosen coefficient of friction. 

With these assumptions the internal friction forces pro- 
duced in the strut depend on the magnitude of the forces on 
the axle, the inclination of the gear, the spacing of the bear- 
ings, and the coefficient of friction between the bearings and 
the cylinder walls. Figure 1(с) schematically illustrates the 
balance of forces acting on the various components of the 
landing gear. The total axial friction in the shock strut is 
the sum of the friction forces contributed by each of the 


. bearings: 


к= (| Fil и ||) 

where | 

Е, axial friction force 

ш coefficient of friction for upper bearing en to 
inner cylinder) 

F, normal force on upper bearing (attached to inner 
cylinder) 

us coefficient of friction for lower bearing (attached to 

outer cylinder) 
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F, normal force on lower bearing (attached to outer 


cylinder) 

[8] factor to indicate sign of friction force 

During the interval prior to the beginning of shock-strut 
motion the friction forces depend on the coefficients of static 
friction; after the strut begins to telescope the coefficients of 
kinetic friction apply. 

From considerations of the balance of moments it can be 
seen from figure 1(c) that 


and ; 
2— 38 
Да 2-12) 
so that 
Fe Ре | Gata) sens | (4) 
J= Е | Na 1 2 1-8 2 
where 
Fy = Fy, sin e— Fy, cos e (4a) 

and з 
Fy, force пога] to strut applied at axle J 
Fy, vertical force applied at axle 

на uorizontal force applied at axle 


ce angle between strut axis and vertical 
l axial distance between upper and lower bearings, for 
fully extended strut 
1. axial distance between axle and lower bearing (attached 
- to outer cylinder), for fully extended strut 

The quantities Fy, Fv, and Fy, are forces applied at the 
axle and differ from the ground reactions by amounts equal 
to the inertia forces corresponding to the respective accelera- 
tion component of the lower mass. Since the inner cylinder 
generally represents only a relatively small fraction of the 
lower mass, the lower mass may reasonably be assumed to be 
-concentrated at the axle. With this assumption, the rela- 
tionships between the forces at the axle and the forces at the 
ground are given by - 


Fa =( Fa és) 


- 


ға 2—1) 


The normal force at the axle can therefore be expressed in 
terms of the ground reactions and the component accelera- 
tions of the lower mass by 


Fy (Fr t aW) sin (2,7 is) cos е (4b) 


where 
Гу, vertical force applied to tire at ground 
Ен, horizontal force applied to tire at ground 
Ws effective mass below shock strut, assumed concentrated 
g at axle 

horizontal acceleration of axle 

vertical acceleration of axle 
In the case of an inclined landing gear having infinite stiff- 
ness in bending, the horizontal displacement of the lower mass 
га1вгеа ед to the vertical displacements of the upper and lower 


Х 
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masses by the kinematic relationship 2,=(2,—2)tan о, as 
previously noted. Double differentiation of this relation- 


ship gives 25— (21-- 24) бап о. Substitution of this expression 


into equation (4b) gives 
: Ws. . ! 
Ен < Еу, sm Ф- Ён, с08 Дата 21 SN e— Ws, sin Ф (4c) 


In equation (4c) the quantity 2, sin e represents the ac- 
celeration of the lower mass normal to the strut axis when the 
gear is rigid in bending. In the case of a gear flexible in 
bending, the normal acceleration of the lower mass is not 
completely determined by the vertical acceleration of the 
upper mass and the angle of inclination of the gear. If it 
should be necessary to take into account, in particular cases, 
the effects of gear flexibility on the relationship between the 
normal force on the axle and the ground reactions, the quan- 
tity 2, sin e in equation (4c) may be replaced by estimated 
values of the actual normal acceleration of the lower mass 88 
determined from consideration of the bending response of the 
gear to the applied forces normal to the gear axis. The effecta 
of gear flexibility are not considered in more detail in the 
present analysis. " 


FORCES ON TIRE 


Figure 2 (а) shows dynamic force-deflection charactoris- 
tics for.a 27-inch smooth-contour (type I) tire inflated to 32 
pounds per square inch. These characteristics were deter- 
mined from time-history measurements of vertical ground 
force and tire deflection in landing-gear drop tests with 8 
nonrotating wheel at several vertical velocities. As can be 
seen, the tire compresses along one curve and unloads along 
another, the hysteresis loop indicating appreciable energy 
dissipation in the tire. There is some question as to whether 


-the amount of hysteresis would be as great if the tire were 


rotating, as in a landing with forward speed. The force- 
deflection curve for a velocity of 11.63 feet per second is for 
& severe impact in which tire bottoming occurs and shows 
the sharp increase in force with deflection subsequent to 
bottoming. 

In figure 2 (b) the same force-deflection characteristics 
are shown plotted on logarithmic coordinates. As can be 
seen, the force exhibits an exponential variation with deflec- 
tion. А systematized representation of the force-deflection 
relationship can therefore be obtained by means of simple 
equations having the form 


T 

Еу == те =m (2) (5) 
where 
Fy, - vertical force, applied to tire at ground 
2 vertical displacement of lower mass from position at 

| initial- contact (radial deflection of tire) 

d overall diameter of tire 
m,r constants corresponding to the various regimes of the 


tire-deflection process 
т’ combined constant, md’ 
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(a) Uniform coordinates. 
Етсовев 2.—Dynamic force-deflection characteristics of tire. 


= 


It may be noted from figure 2 that essentially the same 
force-deflection curve holds during compression for all impact 
velocities, up to the occurrence of tire bottoming, and that in 
figure 2 (b) the slopes -of the curves in each of the several 
regimes of the tire-deflection process are also independent of 
velocity, except in the compression regime following tire 
bottoming. 

Figure 2 also shows simple approximations to the tire 
characteristics which were obtained by fitting straight-line 
segments (long-dashed lines) to the actual force-deflection 
curves in figure 2 (a) for impacts at 8.86 and 11.63 feet per 
second. These approximations, hereinafter referred to as 
linear-segment approximations, are included in a study, 
presented in a subsequent section, to evaluate the degree of 
accuracy required for adequate representation of the tire 
characteristics. The various representations of the tire 
characteristics considered and the pertinent constants for 
each regime of tire deflection are shown in figure 3. 


EQUATIONS OF MOTION 


The internal axial force Fs; produced by the shock strut 
was shown in a previous section to be equal to the sum of the 


+ 


5 : 
Tire deflection, 2, ft. 


(b) Logarithmic coordinates 


hydraulic, pneumatic, and friction forces, as given by 
equation (1). Simce these forces act along the axis of the 
strut, which may be inclined to the vertical by an angle о, 
the vertical component of the axial shock-strut force is given 
by Fs cos e. The vertical component of the force normal 
to the shock strut is given by Fy, sin e. These forces act 
in conjunction with the lift force and weight to produce an 
acceleration of the upper mass. 'The equation of motion 
for the upper mass is 


Fs cos g+ Fy, sin e--L— Wie — 6 (6) 


The vertical components of the axial and normal shock- 
strut forces also act, in conjunction with the weight of the 
lower mass, to produce a deformation of the tire and an 
acceleration of the lower mass. The equation of motion 
for the lower mass is 


Fs cos e+ Fy, sin g+ по а 23— Еу, (22) (7) 


where the vertical ground reaction Ку, is expressed as а. 
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(a) Impact without tire bottoming, Vy,—8.86 feet per second.. 


1 
Tire deflection, 22, ft 


(b) Impact with tire bottoming, Vy,=11.63 feet per second. 


FIGURE 3.—Tire characteristics considered in solutions (logarithmic coordinates). 


function of the tire deflection z}. The relationship between 
Е» : and 2; has been discussed in the previous section on tire 
characteristics. 

By combining equations (0) and (7), the vertical ground 
force can be written in terms of the inertia reactions of the 
upper and lower masses, the lift force, and the total weight. 
The overall dynamic equilibrium is given by 


Ру, (22)= – 그 28947 W (8) 


MOTION PRIOR TO SHOCK-STRUT DEFLECTION 


Conventional oleo-pneumatic shock struts are inflated 
to some finite pressure in the fully extended position. Thus 
the strut does not begin to deflect in an impact until sufficient 
force is developed to overcome the initial preloading imposed 
by the air pressure and internal friction. Since the strut is 


effectively rigid in compression, as well as in bending, 


prior to this instant, the system may be considered to have 
only one degree of freedom during the initial stage of 


the impact. The equations of motion for the one-degreo-of- 
freedom system are derived in order to permit determination 
of the initial conditions required for the analysis of the 
landing-gear behavior subsequent to the beginning of shock- 
strut deflection. 

Since 2,—2,—2 during this first, phase of the impact, 
equation (8) may be written as 


Fr, ()=— 2 $— WK 1) (9) 


where 


Ay 


For the general case of an exponential relationship between 
vertical ground force and tire deflection, equation (5) applies 
and the equation of motion becomes 


И рупа + 7150 (10) 
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The shock strut begins to telescope when the sum of the inertia, weight, and lift forces becomes equal to the vertical 
components of the axial and normal shock-strut forces. At this instant tr, Fs— P; Еу, and equation (6) can be written as 


В Ст) cos ot ly, sin e Ку И —– И} 


— 11) 
Wile . | 
where | | 
Fa, initial air-pressure preload force, ра Аа 
F, static friction at instant t, 
At the instant Е, s=0 and equation (4) becomes 
Еу |Р, Е, (11a) 


where 


к,ан) 2и | 


and ш and из are coefficients of static friction. 


Since the strut is assumed essentially rigid in compression (and also rigid in bending), there is no kinematic displace- 
ment of the lower mass in the horizontal direction up to the beginning of shock-strut deflection, so that 2,==0 and equation 


(4b) becomes 


Wo .. . 
Fy, (Р, РЗ 2,-- Из | sin Ф— Ен, cos o (11b) 
Incorporating equations (Па), (11b), and (9) into equation (11) gives 
Е (ЕК, sin e—cos o) (Kr И — И) – Ен (4K, cos o +sin e) (19) 


г ША 
8 


In equation (12) wherever the + sign appears, the plus signs 
apply when Fy, >0 and the minus signs apply when Fy, <0. 


From equation (10) the vertical displacement of the system 
at the instant t, is given in terms of the corresponding 
acceleration by 


1/7 
„тако 


Integrating equation (10) and noting that 20=0 provides 
the relationship between the vertical velocity and the 
vertical displacement of the system at the beginning of 
shock-strut deflection _ 


jj 22—25 ES of (Е, Ог, | 


In view of the fact that the tre force-deflection curve is 
essentially linear for small deflections, it may be reasonably 
assumed that r--1 for the purpose of determining the time 
after contact at which the strut begins to telescope. With 
this assumption £, can be determined from the relationship 


T Í “42. f j EENET.. LE 
5, d T IE. 27+ тка 1) | 


(13) 


(14) 


(그, sin p— cos о) 


where the general expression for the variable 2 is obtained 
from equation (14) without the subscripts т. Performing 
the indicated integration gives 


t= m [sinc C(1— K;)—sin!GC а-к)": |) (15) 
where | 


C= g 
d ad et — Kg 


The computation of 1, can be greatly simplified by use of 
the following approximation which assumes a linear relation- 
ship between velocity &nd time: 


22, 
i=- 2 
Zot Zr 


(158) 


Equation (158) should be a fairly good approximation in 
view of the relatively short time interval between initial con- 
tact and the beginning of shock-strut motion. 

Equations (12), (13), and (14) permit the determination 
of the vertical acceleration, displacement, and velocity, re- 
spectively, of the system (upper and lower masses) at the 


т = — --. a a MÀ --. 
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beginning of shock-strut deflection. Equation (15) or (15a) permits calculation of the time interval between initial contact 
and this instant. These equations provide the initial conditions required for the analysis of the behavior of the landing 
gear as a system with two degrees of freedom after the shock strut begins to deflect. 

If drag loads are considered, the solution of equation (12) requires knowledge of the horizontal ground force Fy „2t 
the instant t. Since the present analysis does not explicitly treat the determination of drag loads, values of Fyr и have 


to be-estimated, either from other analytical considerations, experimental data, or on the basis of experience. 


MOTION SUBSEQUENT TO BEGINNING OF SHOCK-STRUT DEFLECTION 


Once the sum of the inertia, weight, and lift forces becomes sufficiently large to overcome the preloading force in the 
shock strut due to initial air pressure and internal friction, the shock strut can deflect and the system becomes one having 
two degrees of freedom. Incorporatitig the expressions for the hydraulic, pneumatic, and friction forces (eqs. (2a), (3), 
and (4)) into equation (6) permits the equation of motion for the upper mass to be written as follows: 


Wm. 8 рА ., to п B 1—8 | T 
z арфа S° Paga (Lu) + raj 이 ІС pete | cos y+ K,W— И Fy, sin y=0 (16) 
where | 
Е _ 21-22 
_ совф 
2108 
_ 0080 


and, since Fy = Ру, (22), equation (4c) becomes 


Fy = Еу, (22) sin Ф- Ен, cos er sin o— И sin e 


where Ру (2) is determined {бй the force-deflection characteristics of the tire. 


Fy (а) —mz, as previously noted. 


For the usual type of pneumatic tire, 


Similarly, the equation of motion for the lower mass follows from equation (7): 


E 5 


TA TV стеу rii НЫ [Grad pgs || өвез Б. Дед Раі ө- mc 


(17) 


The overall dynamic equilibrium equation is still, of course, as given by equation (8) 


x 80 в (K,—1)+ Ру, (2:)-0 


Any two of the preceding equations (eqs. (16), (17), and 
(8)) are sufficient to describe the behavior of the landing 
gear subsequent to the beginning of shock-strut motion. 
These equations may be used to calculate the behavior of a 
given landing-gear configuration or to develop orifice and 
metering-pin characteristics required to produce a specified 
behavior for given impact conditions. 'They may also be 
used as & basis for the calculation of dynamic loads in flexible 
airplane structures either by (a) determining the landing- 
gear forcing function-under the assumption that the upper 
mass is a rigid body and then using this forcing function to 
calculate the response of the elastic system or (b) combining 
the preceding equations with the equations representing the 
additional degrees of freedom of the structure; the simul- 
taneous solution of the equations for such a system would 
then take into account the interaction between the deforma- 
tion of the structure and the landing gear. 


SOLUTION OF EQUATIONS OF MOTION 


In the general case the analysis of a landing gear involves 
the solution of the equations of motion given in the section 
entitled “Motion Subsequent to the Beginning of Shock- 
Strut Deflection,” with the initial conditions taken as the 
conditions of motion at the begmning of shock-strut deflec- 
tion, as determined in accordance with the initial impact 
conditions and the equations given in the section entitled 
“Motion Prior to Shock-Strut Deflection.” 


NUMERICAL INTEGRATION PROCEDURES 


In view of the fact that the equations of motion for the 
landing gear subsequent to the beginning of shock-strut 
deflection are highly nonlinear, analytical solution of these 
equations does not appear feasible. In the present report, 
therefore, finite-difference rnethods are resorted to for the 


- step-by-step integration of the equations of motion. Al- 
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though such numerical methods lack the generality of ana- 
lytical solutions and are especially time consuming if the 
calculations are carried out manually, the increasing availa- 
bility of automatic calculating machines largely overcomes 
these objections. 

Most of the solutions presented in this report were obtained 
with а procedure, hereinafter referred to as the “linear pro- 
cedure," which assumes changes in the motion variables to 
be linear over finite time intervals. A few of the solutions 
presented were obtained with a procedure, hereinafter referred 
to as the “quadratic procedure," which assumes a quadratic 
variation of displacement with time for successive intervals. 
The generalized solutions for the simplified equations dis- 
cussed in а subsequent section were obtained by means of 
the Runge-Kutta procedure. The application of these 
procedures is described in detail in appendix А. 


USE OF TIRE FORCE-DEFLECTION CHARACTERISTICS 
In order to obtain solutions for particular cases, it is, of 


course, necessary to have, in addition to information regard- 


ing the physical characteristics of the landing gear, some 
knowledge of the force-deflection characteristics of the tire. 

If extensive data regarding the dynamic tire character- 
istics, such ав shown in figures 2 and 3, are available, an 
accurate solution can be obtained which takes into account 
the various breaks in the force-deflection curves (logarithmic 
coordinates), as well as the effects of hysteresis. In view of 
the fact that the constants m’ and r have the same values 
throughout practically the entire tire compression process 
regardless of the initial impact velocity or the maximum 
load attained, these values of m’ and 7, as determined from 
the force-deflection curves, can be used in the calculation of 
the motion subsequent to the beginning of shock-strut deflec- 
tion until the first break in the force-deflection curve is 
reached prior to the attainment of the maximum force. If 
the conditions for the calculations are the same as those for 
which force-deflection curves are available, the values of 
m’ and r for each of the several regimes subsequent to the 
first break can also be determined directly from the force- 
deflection curves. In general, however, the conditions will 


not be the same and interpolation will be necessary to 


estimate the values of m’ for the subsequent regimes. 
Such interpolation is facilitated, particularly after the maxi- 
mum force-deflection point has been calculated, by the fact 


that each subsequent regime has a fixed value of r, regardless 


of the initial impact conditions. 

The use of the tire-deflection characteristics in the calcula- 
tions is greatly simplified if hysteresis is neglected since the 
values of m’ and r which apply prior to the first break in the 
force-deflection curves are then used throughout the entire 
calculation, except in the case of severe impacts where tire 
bottoming occurs, in which case new values of m/ and r are 
employed in the tire-bottoming regime. А similar situation 
exists with respect to the constants a’ and 0 when the linear 
approximations which neglect hysteresis are used. ‘These 
simplifications would normally be employed when only the 
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tire manufacturer’s static or so-called impact load-deflection 
data are available, as is usually the case. 


EFFECT OF DRAG LOADS 


Although the present analysis permits taking into account 
the effects of wheel spin-up drag loads on the behavior of the 
landing gear, the determination of the drag-load time history 
is not treated explicitly. Thus, if it is desired to consider Ше 
effects of the drag load on the gear behavior, such as in the 
case of a drop test in which drag loads are simulated by 
reverse wheel rotation or in a landing with forward speed, it 
is necessary to estimate the drag load, either by means of 
other analytical considerations or by recourse to experimental] 
data. As a first approximation the instantaneous drag force 
may be assumed to be equal to the vertical ground reaction 
multiplied by a suitable: coefficient of friction и; that is, 
Fg ™= Ру ш, up to the instant when the wheel stops skidding, 
after which the drag force may be assumed equal to zero. 
(The current ground-loads requirements specify a skidding 
coefficient of friction 40.55; limited experimental evidence, 
on the other hand, indicates that и may be as high as 0.7 ог 
as low as 0.4.) In some cases experimental data indicate 
that representation of the drag-load time history can be 
simplified even further by assuming a linear variation of the 
drag force with time during the period of wheel skidding. 

The instant at which the wheel stops skidding can be 
estimated from the simple impulse-momentum relationship 


Io Va, 


rE 


bow 


Fa, а= |” Fy, 04-т------ 


where 
I, polar moment of inertia of wheel assembly about axle 
Рн, initial horizontal velocity 
r4 X radius of deflected tire 
ts, іше of wheel spm-up 
When the drag force is expressed in terms of the vertical 


i 
force, the value of the integral |. Fy, dt can be determined as 


the step-by-step calculations proceed and the drag-force 
term eliminated from the equations of motion after the re- 
quired value .of the integral at the instant of spin-up is 
reached. 


EVALUATION OF ANALYSIS BY COMPARISON OF CALCULATED 
RESULTS WITH EXPERIMENTAL DATA 


In order to evaluate the applicability of the foregoing 
analytical treatment to actual landing gears, tests were 
conducted in the Langley impact basin with a conventional 
oleo-pneumatic landing gear originally designed for a small 
military training airplane. А description of the test specimen | 
and apparatus used is given in appendix В. 

In this section calculated results are compared with ex- 
perimental data for a normal impact and a severe impact 
with tire bottoming. The vertical velocities at the instant 
of ground contact used in the calculations correspond to 
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the vertical velocities measured in the tesis. Equations 
(12), (13), (14), and (15a) were used to calculate the values 
of the variables at the instant of initial shock-strut deflection. 
Numerical integration of equations (16) and (17) provided 
the calculated results for the two-degree-of-freedom system 
subsequent to the beginning of shock-strut deflection. 

In these calculations the discharge coefficient for the orifice 
and the polytropic exponent for the air-compression process 
were assumed to have constant values throughout the impact. 
Consideration of the shape of the orifice and examination 
of data for rounded approach orifices in pipes suggested 
a value of C, equal to 0.9. Evaluation of data for other 
landing gears indicated that the air-compression process 
could be represented fairly well by use of an average value 
of the effective polytropic exponent n=1.12. In view of 
the fact that the landing gear was mounted in а vertical 
position and drag loads were absent in the tests, friction 
forces in the shock strut were assumed to be negligible 
in the calculations. Since the weight was fully balanced 
by lift forces in the tests, the lift factor A, was taken equal to 
1.0. The appropriate exact tire characteristics (See fig. 3) 
were used for each case. | 


NORMAL IMPACT 


Figure 4 presents a comparison of calculated results with - 


experimental data for an impact without tire bottoming at 
a vertical velocity of 8.86 feet per second at the instant of 
ground contact. The exact dynamic force-deflection charac- 
teristics of the tire, including hysteresis, were used in the 
calculations. These tire characteristics are shown by the 
solid lines in figure 2 (a) and values for the tire constants 
m’ and г are given in figure 3 (a). _ 

Calculated time histories of the total force оп the upper 
mass and the acceleration of the lower mass are compared 
with experimental data in figure 4 (а). Similar comparisons 
for the upper-mass displacement, upper-mass velocity, lower- 


` mass displacement, strut stroke, and strut telescoping 


velocity are presented in figure 4 (b). As can be seen, the 
agreement between the calculated and experimental results 
is reasonably good throughout most of the time history. 
Some of the minor discrepancies during the later stages of 
the impact appear to be due to errors in measurement since 
the deviations between the calculated and experimental 
upper-mass accelerations (as represented by the force on 
the upper mass) are incompatible with those for the upper- 
mass displacements, whereas the calculated upper-mass dis- 


placements are necessarily directly compatible with the - 


calculated upper-mass accelerations. The maximum value 
of the experimental acceleration of the lower mass may be 


. somewhat high because of overshoot of the accelerometer. 


In addition to the total force on the upper mass, figure 
4 (a) presents calculated time histories of the hydraulic 
and pneumatic components of the shock-strut force, as 
determined from equations (2) and (3), respectively. It 
can be seen that throughout most of the impact the force 
developed in the shock strut arises primarily from the hy- 
draulic resistance of the orifice. Toward the end of the 
impact, however, because of the decreased - telescoping 
velocities and fairly large strokes which correspond to high 
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compression ratios, the air-pressure force becomes larger 
than the hydraulic force. 


IMPACT WITH TIRE BOTTOMING 


Figure 5 presents a comparison of calculated and experi- 
mental results for a severe impact (Vy,=11.63 ft per sec) 
in which tire bottoming occurred. "The tire force-deflection 
characteristics used in the calculations are shown by the 
solid lines in figure 3 (b). Region (1) of the tire force- 
deflection curve has the same values of the tire constants 
m’ and r as for the case previously discussed. Following 
the occurrence of tire bottoming, however, different values 
of m’ and r apply. These values are given in figure 3 (b). 

It can be seen from figure 5 that the agreement between 
the calculated and experimental results for this case is 
similar to that for the comparison previously presented. 

The calculated instant of tire bottoming is indicated in 
figure 5. When tire bottoming occurs, the greatly increased 
stiffness of the tire causes a marked increase in the shock- 
strut telescoping velocity, as is shown in the right-hand 
portion of figure 5 (b). Since the strut is suddenly forced to 
absorb energy at a much higher rate, an abrupt increase 
in the hydraulic resistance takes place. The further increase 
in shock-strut force immediately following the occurrence of 
tire bottoming 1з evident from the left-hand portion of 
figure 5 (a). The sudden increase in lower-mass acceleration 
at the instant of tire bottoming can also be seen. 

In this severe impact the hydraulic resistance of the orifice 
represents an even greater proportion of the total shock- 
strut force than was indicated by the calculated results for 
an initial vertical velocity of 8.86 feet per second previously 
discussed. | i 

The foregoing comparisons indicate that the analytical 
treatment presented, in conjunction with reasonably straight- 
forward assumptions regarding the parameters involved in 
the equations, provides a fairly accurate representation of 
the behavior of a conventional oleo-pneumatic landing gear. 


PARAMETER STUDIES 


In the previous section comparisons of calculated results 
with experimental data showed that the equations which’ 
have been developed provide a fairly good representation of 
the behavior of the landing gear for the impact conditions 
considered. In view of the fact that the equations are 
somewhat complicated and require numerical values for 
several parameters such as the tire force-deflection constants 
т’ and т, the orifice discharge coefficient Си, and Ше poly- 
tropic exponent n, which may not be readily or accurately 
known in the case of practical engineering problems, it 
appears desirable (a) to determine the relative accuracy 
with which these various parameters have to be known and 
(b) to investigate the extent to which the equations can be 
simplified and still yield useful results. In order to accom- 
plish these objectives, calculations have been made to 
evaluate the effect of simplifying the force-deflection charac- 
teristics of the tire, as well as to determine the effects which 
different values of the orifice discharge coefficient and the 
effective polytropic exponent have on the calculated behavior. 
The results of these calculations are discussed in the present 
section. The question of simplification of the equations of 
motion is considered in more detail in a subsequent section. 
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Figure 4.—Comparisons between calculated results and experimental data for normal impact; 
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Figure 5.— Comparisons between calculated results and experimental data for impact with tire bottoming; solution with 
exact exponential tire characteristics. Vy =11.63 feet per second; C4—0.9; n—1.12. | 
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(b) Time histories of landing-gear velocities and displacements. 
Figure 5.—Concluded. 
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REPRESENTATION OF TIRE FORCE-DEFLECTION CHARACTERISTICS 


In order to evaluate the degree of accuracy required for 
adequate representation of the tire forcé-deflection charac- 
teristics, comparisons are made of the calculated behavior 
of the landing gear for normal impacts and impacts with 
tire bottoming when the tire characteristics are represented 
in various ways. First, the force-deflection characteristics 
wil be assumed to be exactly as shown by the solid-line 
curves in figure 2 (b), including the various breaks in the 
curve and the effects of hysteresis. These characteristics 
are referred to hereinafter as the exact exponential tire 
characteristics. The effects of simplifying the representa- 
tion of the tire characteristics will then be investigated by 
considering (a) the exponential characteristics without 
hysteresis; that is, the tire will be assumed to deflect and 
unload along the same exponential curve, (b) the linear- 
segment approximations to the tire characteristics (long- 
dashed lines), which also neglect hysteresis, and (c) errors 
introduced by neglecting the effects of tire bottoming m the 
case of severe impacts. The calculated results presented in 
this study make use of the relationships between vertical 
force on the tire and tire deflection, as shown in figures 
3 (a) and 3 (b). | 

Figure 6 presents a comparison of the calculated results 
for а normal impact at a vertical velocity of 8.86 feet рег 
second, whereas figure 7 permits comparison of the solutions 
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for a severe impact, involving tire bottoming, at a vertical 
velocity of 11.63 feet per second. In figures 6 and 7 the 
solid-line curves represent solutions of the landing-gear 
equations when the exact exponential relationships between 
force and tire deflection are considered. Since these solu- 
tions were previously shown to be in fairly good agreement 
with experimental data (figs. 4 and 5), they are used as a 


basis for evaluating the results obtained when tire hysteresis 


is neglected and the force-deflection characteristics are repre- 
sented by either simplified exponential or linear-segment 
relationships. 

As in the calculations previously described, the solutions 
were obtained in two parts. During the first stage of the 
impact the shock strut was considered to be rigid until 
sufficient force was developed to Overcome the initial air- 
pressure force. The calculations for the landing-gear behav- 
ior subsequent to this instant were based on the equations 
which consider the gear to have two degrees of freedom. 
Time histories of the upper-mass acceleration calculated on 
the basis of a rigid shock strut are shown by the dotted 
curves in figures 6 and 7. ‘These solutions show the greatest 
rate of increase of upper-mass acceleration possible with 
the exponential tire force-deflection characteristics con- 
sidered. Comparison of these solutions with those for the 
two-degree-of-freedom system indicates the effect -of the 
shock strut in attenuating the severity of the impact. . 
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(a) Time histories of upper-mass acceleration and lower-mass acceleration. 


FigvRE 6 — Effect of tire characteristics on calculated landing-gear behavior in normal impact. 
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` (b) Time histories of landing-gear displacements and velocities. 
FIGURE 6.— Continued. | 
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(c) Time histories of shock-strut stroke and velocity. 
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Етаове 7.—LEffect of tire characteristics on calculated landing-gear behavior for impact with tire bottoming. 
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Vy, 7711.63 feet per second; Cz=0.9; n—1.12. 
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(b) Time histories of landing-gear displacements and velocities 
Етсовв 7.—Continued. 
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(c) Time histories of shock-strut stroke and velocity. А 
FIGURE 7.—Concluded. 


Normal impact.—In the case of the normal impact at а 
vertical velocity of 8.86 feet per second, figure 6 shows that 
the solution obtained with the exponential force-deflection 
variation which neglects hysteresis and the solution with the 
linear-segment approximation to the tire characteristics are 
in fairly good agreement with the results of the calculation 
based on the exponential representation of the exact tire 
characteristics. The greatest differences between the solu- 
tions are evident in the time histories of upper-mass and 
lower-mass accelerations; considerably smaller differences are 
obtained for the lower-order derivatives, as might be ex- 
pected. With regard to the upper-mass acceleration, the 
- three solutions are in very good agreement during the early 
stages of the impact. In the case of the simplified exponen- 
tial characteristics, neglect of the decreased slope of the. 
force-deflection curve between the first break and the maxi- 
mum (regime (2) in fig. 3 (2)) resulted in the calculation of a 
somewhat higher value of the maximum upper-mass acceler- 
ation than was obtained with the exact tire characteristics. 
For the simplified exponential and linear-segment character- 
istics, neglect of hysteresis resulted in the calculation. of 
somewhat excessive values of upper-mass acceleration sub- 
sequent to the attainment of the maximum vertical load. 
It is of interest to note that the calculated results for the 
exponential and linear-segment characteristics without hys- 
teresis were generally in quite good agreement with each 
other throughout the entire duration of the impact, although 
the assumption of lmear-segment tire force-deflection char- 
acteristics did result in somewhat excessive values for the 
maximum lower-mass acceleration. On the whole, . the 
simplified tire force-deflection characteristics considered per- 
mit calculated results to be obtained which represent the 
behavior of the landing gear in normal impacts fairly well. 


Impact with tire bottoming.—In the case of the severe 
impact at a vertical velocity of 11.63 feet per second, the 
effects of tire bottoming on the upper-mass acceleration, the 


 lower-mass acceleration, and the strut telescoping velocity 


are clearly indicated in figure 7 by the calculated results 
based on the exact tire characteristics. As can be воеп, the 
linear-segment approximation to the tire deflection character- 
istics which takes into account the effects of tire bottoming 
resulted in a reasonably good representation of the landing- 
gear behavior throughout most of the time history. On the 
other hand, as might be expected, the calculations which 
neglected the effects of bottoming on the tire force-deflection 
characteristics did not reveal the marked increase in the 
upper-mass acceleration due to the increased stiffness of tho 
tire subsequent to the occurrence of bottoming. It is also 
noted that the discrepancies in the calculated upper-mass 
acceleration due to neglect of hysteresis in the later stages 
of the impact are more pronounced in this case than in the 
impact without tire bottoming previously considered, again 
as might be expected. 


EFFECT OF ORIFICE DISCHARGE COEFFICIENT 


In view of the fact that there is very little information 
available regarding the magnitude of discharge coefficients 
for orifices in landing gears, it appears desirable to evaluate 
the effect which differences in the magnitude of the orifice 
coefficient can have on the calculated results. Figure 8 
presents comparisons of calculated results for а range of 
values of the orifice discharge coefficient C, between 1.0 
and 0.7. Тһе four solutions presented are for Ше same set 
of initial conditions as the normal impact without tire 
bottoming previously considered and are based on the 
exponential tire force-deflection characteristics which neglect 


hysteresis. 
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Fiaurs 8.—Effect of orifice discharge coefficient on landing-gear behavior; calculations with exponential tire characteristics 
without hysteresis. Vy,=8.86 feet per second; n=1.12. 
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(b) Time histories of landing-gear displacements and velocities, 
Етаунв 8.—Continued. 
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Етатувв 8.—Concluded. 


These calculations show that a decrease in the orifice dis- 
charge coefficient results in an approximately proportional 
increase in the  upper-mass acceleration. This vari- 
ation is to be expected since the smaller coefficients cor- 
respond to reduced effective orifice areas which result in 
greater shock-strut forces due to increased hydraulic 
resistance. As a result of the increased shock-strut force 
acting downward on the lower mass, the maximum upward 
acceleration of the lower mass is reduced with decreasing 
values of the discharge coefficient. The increase in shock- 
strut force with decreasing discharge coefficient also results 
in a decrease in the strut stroke and telescoping velocity but 
an increase in the lower-mass velocity and displacement, 
as might be expected. However, since the increases in 
lower-mass displacement and velocity are smaller than the 
decreases in strut stroke and telescoping velocity, the upper- 
mass displacement and velocity are reduced with decreasing 
orifice discharge coefficient. 

These comparisons show that the magnitude of the orifice 
coefficient has an important effect on the behavior of the 
landing gear and indicates that a fairly accurate determi- 
nation of the numerical value of this parameter is necessary 
to obtain good results. Я 


EFFECT ОР AIR-COMPRESSION PROCESS 
Since the nature of the air-compression process in a shock 
strut is not well-defined and different investigators have 


assumed values for the polytropic exponent ranging any- 
where between the extremes of 1.4 (adiabatic) and 1.0 


(isothermal), it appeared desirable to evaluate the im- 
portance of the air-compression process and to determine the 
extent to which different values of the polytropic exponent 
can influence the calculated results. Consequently, solutions 
have been obtained for three different values of the poly- 
tropic exponent; namely, n=1.3, 1.12, and 0. 

The value 1— 1.3 corresponds to a very rapid compression 
in which an adiabatic process is almost attained. ‘The 
value n=1.12 corresponds to a relatively slow compression 
in which the process is virtually isothermal. The value 
7,—0 is completely fictitious since it implies constant air 
pressure within the strut throughout the impact. The 
assumption n=0 has been considered since it makes one of 
the terms in the equations of motion & constant and permits 
simplification of the calculations. The three solutions 
presented are for the same set of initial conditions as the 
normal impact without tire bottoming previously con- 
sidered and are based on the exponential tire force-deflection 
characteristics which neglect hysteresis. 

Figure 9 shows that the air pressure contributes only г, 
relatively small portion of the total shock-strut force through- 
out most of the impact since the compression ratio is rela- 
tively small until the later stages of the impact. ‘Toward 
the end of the impact, however, the air-pressure force 
becomes a large part of the total force since the compression 
ratio becomes large, whereas the hydraulic resistance de- 
creases rapidly as the strut telescoping velocity is reduced 
to zero. 
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(a) Time histories of upper-mass force and lower-mass acceleration. 


Гісунк 9,—Effect of polytropic exponent; calculations with exponential tire characteristics without hysteresis, 
Vr =8.86 feet per second; Cz=0.9. 
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(b) Time histories of landing-gear displacements and velocities. 
Fiaurn 9.—Continued. 
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(c) Time histories of shock-strut stroke and velocity. 
Етсовв 9.—Concluded. 


As a result, the calculations show that the magnitude of 
the polytropic exponent has only a very small effect on the 
behavior of the landing gear throughout most of the impact. 
For the practical range of polytropic exponents, variations 
in the air-compression process result in only minor differ- 
ences in landing-gear behavior, even during the very latest 
stages of the impact. The assumption of constant air 
pressure in the strut throughout the impact (п-<0), however, 
does lead to the calculation of excessive values of stroke and 
of the time to reach the maximum stroke. The time history 
of the shock-strut force calculated on the basis of this 
assumption is, on the other hand, in quite good agreement 
with the results for the practical range of air-compression 
processes. 

On the whole it appears that the behavior of the landing 
gear is relatively insensitive to variations in the air- 
compression process. The foregoing results suggest that, in 
many cases, fairly reasonable approximations for the landing- 
gear force-time variation might be obtained even if the air- 
pressure term in the equations of motion were completely 
neglected. 


SIMPLIFICATION OF EQUATIONS OF MOTION 


The preceding studies have indicated that variations in the 
tire force-deflection characteristics and in the air-compression 
process individually have only a relatively minor effect on the 
calculated behavior of the landing gear. These results sug- 
gest that the equations of motion for the landing gear might 


be simplified by completely neglecting the internal air- | 


pressure forces in the shock strut and by considering the tire 
force-deflection characteristics to be linear. With these as- 
sumptions, the equations of motion for the upper mass, 
lower mass, and complete system (eqs. (16), (17), and (8)) can 
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be written as follows for the case where the wing lift is equal 
to the weight and the internal friction is neglected: 


= Mi atA 20 Ру 0 
Wa 5 
F3 28 A (21— Z3) +az:+b— W:=0 (1 8) 
=" ын Satazatb= 0 
where | 
8 ра, 
ЗІН 2 (СА „Жсовф 
and 
а slope of linear approximation to tire force-deflection 
characteristics 


6 value of force corresponding to zero tire deflection, as 
determined from the linear-segment approximation to 
the tire force-deflection characteristics 

The motion variables at the beginning of shock-strut do- 
flection can be readily determined in a manner similar to 
that employed in the more general treatment previously dis- 
cussed. For the simplified equations the variables at the 
instant t, are given by 


(19) 


In most cases the term T z,* is small in comparison with 242 


so that 2, 20. 
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(a) Time histories of upper-mags acceleration and lower-mass acceleration. 
Fiaury 10.—Evaluation of calculated results for simplified systems. Vy,=8.86 feet per second; Са==0.9. 


The values determined from equations (19) аге used as 
initial conditions in the solution of equations (18). 

The fact that the lower mass is a relatively small fraction of 
the total mass suggests that the system might be simplified 
even further without greatly modifying the calculated results 
by assuming the lower mass to be equal to zero. With this 
assumption f,—0 and the initial values of the variables in 
equations (18) correspond to the conditions at initial contact. 


EVALUATION OF SIMPLIFICATION 


In order to evaluate the applicability of these simplifica- 


tions, the behavior of the landing gear has been calculated in | 


accordance with equations (18) for an impact with an initial 
vertical velocity of 8.86 feet per second. А similar calcula- 
tion has been made with the assumption ]/;—0. These 
results are compared in figure 10 with the more exact solu- 
tions previously presented in figure 4, which include con- 
sideration of the air-compression springing and the exact 
exponential tire characteristics. A time history of the lower- 
mass acceleration is not presented for the case where Из is 
assumed equal to zero since the values of 2:/g have no 
significance in this case. | 

Figure 10 shows that the two simplified solutions are in 
quite good agreement with each other, as might be expected, 
and are also in fairly good agreement with the more exact 
results. Neglecting the air-pressure forces and assuming а 
linear tire force-deflection variation resulted in the calcula- 
tion of slightly lower values for the maximum upper-mass 
acceleration and somewhat higher values for the maximum 
stroke than were obtained with the more exact equations. 


The effect of neglecting the lower mass was primarily to 
reduce the lower-mass displacement (tire deflection), as a 
result of the elimination of the lower-mass inertia reaction. 
On the whole, it appears that the assumptions considered 
permit appreciable simplification of the equations of motion 
without greatly impairing the validity of the calculated results. 


GENERALIZED TREATMENT 

Equations and solutions.—By writing the simplified equa- 
tious of motion in terms of dimensiouless variables, general- 
ized solutions can be obtained for a wide range of landing-gear 
and impact parameters which may be useful in pre- 
liminary design. If W, is taken equal to zero and it is 
further assumed that the tire force-deflection curve is 
represented by а single straight line through the origin 
(b—0 throughout the impact), equations (18) reduce to 


m #1-Е А(21— 23) —0 


A (21— as —a23—0 (20) 


- 21-ага:- 


+ 


where И A, апа а аге constants, as previously defined, and 


any two of the foregoing equations are sufficient to describe 
completely the behaviorof thesystem. With this representa- 
tion of the system, the shock strut begins to deflect at the in- 
stant of initial contact (¢,=0). Thus, the initial conditions 
for equations (20) are the initial impact conditions; namely, 
21) 23,70 and 21, 23,77 2. 
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(c) Time histories of shock-strut stroke and velocity. 
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ANALYSIS OF LANDING-GEAR BEHAVIOR 


As can be seen, with the equations in this form, the solution 
depends on five parameters, namely, m, A, a, and the 


initial conditions 26 and 2. However, since 2,—0 in all 
cases, the number of variable parameters is reduced to four. 
In view of the fact that these parameters are independent of 
one another and each may take оп a large range of values, a 
great many solutions and a large number of graphs would be 
required to cover the entire range of landing-gear and impact 
parameters with the equations in the form of equations (20). 

The number of independent parameters which have to be 
considered may be greatly decreased by the introduction of 
generalized dimensionless variables and the corresponding 
transformations of equations (20). In this case, generalized 
variables can be obtained which permit transformation of the 
equations of motion to a form which does not rovolve any 
constants. With the equations iu this form, there is only 
one variable parameter, namely, the initial velocity paramı- 
eter. To determine the generalized variables which satisfy 
the aforementioned requirements, let 


U=Za 
and 
0--%8 
Thus, . 
и" Hu. z= 
0 В 
and 


Substituting these new variables permits equations (20) 
to be written as 


аға а И то. 
Uy p ы” =0 


(20а) 


pec 


The number of independent parameters will be reduced 
if all the combined constants in equations (20a) are set 
equal to one another, that is, let 


а[8* 
— Wilg 


From this relationship, it can be seen that 


ac A 
Wig 


| a 
В Тр; 


Thus the generalized variables become 


(uu) ә rg) 


аа а 


321605—50——70 


and 


1105 
and 


-a (4) 
WW 


With these new variables equations (20) can be written as? 


where 


(а — 9)? —0 


(uy! — 4)? —u5—0 (21) 


uy’ +u=0 


where any two of these equations are sufficient to describe 
the behavior of the system. | | 

Inasmuch ав equations (21) do not involve any constants, 
their solutions are completely determined by the initial val- 
ues of the variables. Since the displacements at initial 
contact W, and w are equal to zero and the initial velocities 
Шо and u’ are equal, the only parameter is the initial 
dimensionless velocity 


"EPI ES: 
и А! Wia 


where tp’ =, —us,'. 


Generalized solutions of equations (21) are presented in 
figure 11 for values of 10 corresponding іо a wide range of 


 landing-gear and impact parameters. Parts (а) to (e) of 


figure 11 show the variations of the dimensionless variables 
during the impact; parts (f) and (g) show the maximum 
values of the more important variables as functions of 2)’. 
Part (h) shows the shock-strut effectiveness з, and the 
landing-gear effectiveness пу. The shock-strut effective- 
ness, sometimes called ''efficiency" and, in Europe, ''plani- 
metric ratio," is defined as 


where о--шщ--шз is the dimensionless shock-strut stroke, 
Since 7, represents the ratio of the energy actually absorbed 
by the shock strut to the maximum energy which the strut 
could possibly absorb for any combination of maximum 
acceleration (or load) and maximum stroke, 15 serves as a 
measure of the extent to which a given combination of 
maximum load and stroke has been utilized to absorb the 
energy of an impact. А similar measure of the energy 
absorption effectiveness of the landing gear as a whole is 
given by т: which is defined by 
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Я 44” au, 
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3 Equations (21) may be reduced to a single equation in one variable by differentiating the 
last equation and substituting for uy in the firstequation. This gives (u+ Hu” m0, 

By Introducing the new variable юни”, this equation may be reduced to the second- 
order equation (to-]-15/)1-]-t5' 0, subject to the Initial conditions козиф and ку 20, 
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(a) Relationship between upper-mass acceleration, lower-mass displacement, and time. 
Fraunz 11.—Generalized solutions for simplified system, 
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(b) Relationship between upper-mass displacement and time. 
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Dimensionless time, 0 
(c) Relationship between upper-mass velocity, lower-mass velocity, and time. 
Fiaurs 11.—Continued. 


The generalized results presented in figure 113 can be 
used to estimate the performance of a given landing gear of 
known configuration for particular impact conditions or to 
choose the dimensions for a lending gear when the impact 
conditions and desired performance are specified. 

Applicability of solutions.—To illustrate the applicability 
of the generalized solutions, the curves of figure 11 have 
been applied to the previously considered case of the normal 
impact at an initial vertical velocity of 8.86 feet per second 
for comparison with the more exact solution presented in 
figure 4. In order to make use of the generalized solutions 

3 Although time-history solutions are presented for values of uo’ as small as 0.5, it will be 
noted that values of ti mapt Tmas, з, and y,, are not given for values of tt «1.5. It can bo 
seen from the time histories that the characteristics of the solution in the later stages of the im- 
pact change as uo’ becomes small; іп particular, t; increases and tle curve of timas ав 8 func 
tion of to’ appears to reach a minimum at some value of 1051.5. Furthermore, the later 
Stages of the solutions greatly stretch out in time and appear to be almost asymptotic in 
character, Several different analytical, numerical, and analogue methods were applied In an 


attempt to study this phase of the problem further but the extremely slow rate of change of 
the variables in this region prevented successiul completion of Ше solutions. 


it is first necessary to approximate the tire force-deflection 
characteristics by a simple linear variation. Two such 
linear approximations which might be considered suitable 
for this purpose are shown in figure 12. Linear approxi- 
mation І is a straight line through the origin having a slope 
a=18.5X 10° pounds per foot (@’=ad=41.6X10? Ib). This 
value of а and the other pertinent landing-gear and impact 
parameters result in a value of the initial dimensionless 
velocity parameter u,’ =2.57. Linear approximation П is 
a straight line with slope а--21.8х108 pounds per foot 
(a’=47.910? Ib) which does not pass through the origin 
but intersects the displacement axis at a value of 
Zary су 0-0508 foot. With this value of а, му —2.39. 


Since the solutions of figure 11 have been calculated only 
for integral values of ug’, curves for the foregoing values of 
"uy were graphically interpolated by cross plotting. These 
results were then converted to dimensional values by multi- 
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(d) Relationship between upper-mass acceleration and strut stroke. 
Ғүсовв 11.— Continued. 


plying the dimensionless variables by the appropriate 
constants. ‘The results obtained are compared in figure 13 
with the more exact solution presented in figure 4. The 


values based on linear approximation I have been plotted © 


exactly as determined from the generalized solutions. The 


results for linear approximation П, however, have been’ 


displaced relative to the origin of coordinates as indicated 

in the following discussion. | 

The assumption of linear approximation I implies that 

the system must move a distance equal to ау ut) after 
Е 


initial contact (at constant velocity since the wing lift is 
taken equal to the weight) before any finite ground reaction 
can develop. The derivation of the equations of motion, 
on the other hand, assumes that the ground reaction in- 
creases linearly with deflection from the instant of initial 
contact. As a result, the equations of motion do not apply 


until after the system has attained а displacement equal to ' 


23 ғу,-оу which occurs at a time after initial contact 
t 


22 " 
t= Cre 2 In other words, the equations of motion 
V 


0 
apply to a coordinate system transformed so that the tire 


force-deflection relationship passes through the origin; that 
is, & coordinate system displaced by бк, -6 relative to 
9. 


Led 


the coordinate system originating at the point of initial 

contact. It therefore follows that the upper-mass and 

lower-mass displacements determined from the generalized 

solutions for the case of linear approximation II must be 

increased by & constant amount equal to 2 ғу „o. IN this 
£ 


case 0.0508 foot, and all results must be displaced in 
ғу) 
y 


; in this caso 


time by a constant increment At= 
Vo 


2-0508 0.0057 second, relative to the instant of initial 
contact. These corrections have been incorporated in plot- 


ting the curves for linear approximation II shown in figure 13. 
As can be seen, the results obtained by application of the 
generalized solutions, particularly by the method employing 
linear approximation 1I, are in fairly good agreement with the 
more exact solution. The discrepancies which exist are 
attributable to the neglect of the shock-strut preloading and 
springing provided by the air-pressure force, neglect of the 
lower mass, and to differences between the very simple tire 
force-deflection relationships assumed and the exact tire 
characteristics. On the whole, it appears that the general- 
ized results offer a means for rapidly estimating the behavior 
of the landing gear within reasonable limits of accuracy and 
may therefore be useful for preliminary design purposes. 
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(e) Relationship between upper-mass acceleration and upper-mass прасе 
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SUMMARY OF RESULTS AND CONCLUSIONS 


A theoretical study has been made of the behavior of the 
conventional type of oleo-pneumatic landing gear during the 
process of Janding impact. The basic analysis is presented 
in a general form and treats the motions of the landing gear 
prior to and subsequent to the beginning of shock-strut 
deflection. In the first phase of the impact the landing gear 
is treated as a single-degree-of-freedom system in order to 
determine the conditions of motion at the instant of initial 
shock-strut deflection, after which instant the landing gear is 
considered as a system with two degrees of freedom. ‘The 
equations for the two-degree-of-freedom system consider 
such factors as the hydraulic (velocity square) resistance of 
the orifice, the forces due to air compression and internal 
friction in the shock strut, the nonlinear force-deflection 
characteristics of the tire, the wing lift, the inclination of the 
landing gear, and the effects of wheel spin-up drag loads. 

The applicability of the analysis to actual landing gears 
has been investigated for the particular case of a vertical 
landing gear in the absence of drag loads by comparing 
calculated results with experimental drop-test data for corre- 
sponding impact conditions, for both a normal impact and a 
severe impact involving tire bottoming. 

Studies have also been made to determine the effects of 
variations in such parameters as the dynamic force-deflection 


characteristics of the tire, the orifice discharge coefficient, and 
the effective polytropic exponent for the air-compression 
process, which might not be known accurately in practical 
design problems. 

In addition to the more exact treatment an investigation 
has also been made to determine the extent to which the 
basic equations of motion can be simplified and still yield 
useful results. Generalized solutions of the simplified 
equations obtained are presented for a wide range of landing- 
gear and impact parameters. 

On the basis of the foregoing studies the following con- 
clusions are‘indicated: 

1. The behavior of the landing gear as calculated from the 
basic equations of motion was found to be in good agreement 
with experimental drop-test data for the case of a vertical 
landing gear in the absence of drag loads, for both a normal 
impact and a severe impact involving tire bottoming. 

2. Astudy of the effects of variations in the force-deflection 
characteristics of the tire indicates that 

a. In the case of & normal impact without tire bottoming, 
reasonable variations in the force-deflection characteristics 
of the tire have only a relatively small effect on the calculated 
behavior of the landing gear. Approximating the rather 
complicated force-deflection characteristics of the actual tire 
by simplified exponential or linear-segment variations appears 
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(f) Variation of maximum upper-mass acceleration and time to reach maximum upper-mass acceleration with initial velocity parameter. 
| Етсовв 11.—Continued. à 


to be adequate for practical purposes. Tire hysteresis was 
found to be relatively unimportant. 

b. In the case of a severe impact involving tire bottoming, 
the use of simplified exponential and linear-segment approxi- 
mations to the actual tire force-deflection characteristics 
which neglect the effects of tire bottoming, although adequate 
up to the instant of bottoming, fails to indicate the pro- 
nounced increase in landing-gear load which results from 
bottoming of the tire. The use of exponential or linear- 
segment approximations to the tire characteristics which 
take into account the increased stiffness of the tire that re- 
sults from bottoming, however, yields good results. 

3. A study of the importance of the discharge coefficient 
of the orifice indicates that the magnitude of the discharge 
coefficient has a marked effect on the calculated behavior of 
the landing gear; a decrease in the discharge coefficient (or 
the product of the discharge coefficient and the net orifice 
area) résults in an approximately proportional increase in 
the maximum upper-mass acceleration. 

4. A study of the importance of the air-compression process 
in the shock strut indicates that the air springing is of only 
minor significance throughout most of the impact, and that 
variations in the effective polytropic exponént n between the 
isothermal ‘value of 1.0 and the near-adiabatic value of 1.3 
have only a secondary effect on the calculated behavior of 


the landing gear. Even the assumption of constant air 
pressure in the strut equal to the initial pressure (n=0) 
yields fairly good results, which may be adequate for many 
practical purposes. 

5. An investigation of the extent to which the equations of 
motion for the landing gear can be simplified and still yield 
acceptable calculated results indicates that, for many prac- 
tical purposes, the air-pressure force in the shock strut can 
be completely neglected, the tire force-deflection relationship 
can be assumed to be linear, and the lower or unsprung mass 
can be taken equal to zero. 

6. Generalization of the equations of motion for the 
simplified system described in the preceding paragraph 
shows that the behavior of this system is completely deter- 
mined by the magnitude of one parameter, namely, the 
dimensionless initial-velocity parameter. Solution of these 
generalized equations in terms of dimensionless variables 
permits compact representation of the behavior of the system 


. for a wide range of landing-gear and impact parameters, 


which may be useful for rapidly estimating landing-gear 
performance in preliminary design. 


LANGLEY ABRONAUTICAL LABORATORY, 
NATIONAL ADVISORY CoMMITTEE FOR ÁERONAUTICS, 
LawaLEx Етегі, Va., Мау 1, 1952. 
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(g) Variation of maximum upper-mass displacement, maximum lower-mass displacement, and maximum strut stroke with initial velocity parameter. 
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APPENDIX A 
NUMERICAL INTEGRATION PROCEDURES 


As previously noted, most of the specific solutions presented 
in this report were obtained with a numerical integration 
procedure, termed the ‘linear procedure," which assumes 
changes in the variables to be linear over finite time intervals. 
With this procedure a time interval e—0.001 second was 
used in order to obtain the desired accuracy for the particular 
cases considered. A few of the specific solutions presented 
were obtained by means of a procedure, termed the “quad- 
ratic procedure,” which assumes a quadratic variation of 
displacement with time for successive intervals. This pro- 
cedure, although requiring somewhat more computing time 
per interval, may permit an increase in the interval size for 
а given accuracy, in some cases allowing a reduction in the 
total computing time required. In the case of the more 
exact equations of motion the accuracy of the quadratic pro- 
cedure with a time interval of 0.002 second appears to be 
equal to that of the linear procedure with an interval of 
0.001 second. Although the accuracy naturally decreases 
with increasing interval size, the loss in accuracy for pro- 
portionate increases in interval size appears to be smaller for 
the quadratic than for the linear procedure. In the case of 
the simplified equations of motion reasonably satisfactory 
results were obtained in test computations with the quadratic 
procedure for intervals as large as 0.01 second, whereas the 
linear procedure was considered questionable for intervals 
larger than 0.002 second. І 

The generalized solutions presented, because of Ње 
relatively simple form of the equations of motion, were 
obtained with the well-known Runge-Kutta procedure. 
A study of the allowable interval size resulted in the use of 
an interval Аб--0.08, which corresponds to а time interval of 
about 0.005 second for the landing gear under consideration. 

LINEAR PROCEDURE 

In this step-by-step procedure the variations in dis- 
placement, velocity, and acceleration are assumed to be 
linear over each finite time interval e. The method, as 
used, involves one stage of iteration. Linear extrapolation 
of the velocity at the end of any interval is used to obtain 


estimated values of velocity and displacement for the next 


interval. These values are then used to calculate values of 
the acceleration in accordance with the equations of motion. 
Integration of the acceleration provides improved values of 
the velocity and, if desired, the displacement and accelera- 
tion. In this procedure all integrations are performed by 
application of the trapezoidal rule. 

The following derivation illustrates the application of the 
linear procedure to the equations of motion for the landing 
gear, which apply subsequent to the beginning of shock-strut 
deflection at time і. Іп the example presented internal 
friction forces are neglected in order to simplify the deriva- 

1114 


. tion. 


However, the same general procedure can be used if 
these, or other complicating effects, are included in the 
equations. 

For the case under consideration the equations of motion 
(eqs. (16), (17), and (8)) can be written as follows: 


DÀ atA В —Oler— 29)" +D=0 (Al) 


u 22— А(21— 22. — B[1 — 0(21— 23)] ^" + Fy (23) — Wa=0 
(A2) 
P a pP 5-Р, (e+ E=0 (A3) 
where . 
^ 
[3] 2(G,A,)*cose 
В--ра Аасовф 
_ 4 
00С08ф 
D=K, W-— WA 
“ E= W(Kr— 1) 
Solving equation (АЗ) for 21 gives 
i= [F—@3,—HFy, (23) r | (A4) 
where 
Ww 
—W, (1— Ky)g 
- 
- 


Integrating equation (A4) with respect to ¢ between the 
limits 4, and ¢ and noting that 2, = =z, gives 
КМК ОЕ. кыы, | "Fy (ейт (A5) 


where т=(1—#,). 
Integrating again and noting that 2; <<2а << г, gives 


| | | 때 
zy (144) ва -бга--Н |, |, Fy (ет dr 
E (А6) 
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Substituting for 2; and 2; in equation (A3) gives 


“ы (4 (а S yocp ЕО ен |, ржораи 


в1-о|а+фечзи-29+ РН | | Fre (0007 dr |} "р, 2+) | 


The motion of the landing gear subsequent to the beginning 
of shock-strut deflection is determined by means of a step- 
by-step solution of equation (A7). This numerical procedure 
yields time histories of the lower-mass motion variables 
г, 2, and 2, from which the motion variables for the upper 
mass 21, 21, and 2, can be calculated by means of equations 
(A4), (A5), and (A6). 

The initial conditions for the step-by-step procedure are 


21,20 22620 ^, 
2000 228.0 Žr (A8) 
214-0 2, 2977 Er 


where г, 2, and 2, are the conditions of motion at the 
beginning of shock-strut deflection as determined from the 
solution for the one-degree-of-freedom system. 

Estimated values of the lower-mass velocity at the end of 
the first time increment e following the beginning of shock- 
strut deflection can be obtained from the expression 


2 per tees (A9) 
or, a3 a first approximation, 
22,411 2: 
The corresponding displacement is given by 
2 1 =2,+5 (22,4 2.) (А10) 


After the initial conditions and the conditions at the end 
of the first time increment are established, a step-by-step 
calculation of the motion can be obtained by routine opera- 
tions as indicated by the following general procedure which 
applies at any time r=ne after the beginning of the process. 
The operations indicated are based on integration by appli- 
cation of the trapezoidal rule: 

2 +5 (2а, 12а, ) (АП) 


= 2231102231 2233) = Za, 1 


2 
21, == га 1 На (22.1 Е 22 p= Zan T€ had (22, t Za, 3) 
(A12) 


With the estimated values # and zf, the acceleration of 
the lower mass can be determined by substitution in the 


A ul 
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(A7) 


appropriate integrodifferential equation for the system, 
equation (A7) in the present case. Thus 
Za =J (22, 22, та) (A13) 
Ín equation (À7) the integral expressions can also be evalu- 
ated by application of the trapezoidal rule. For example, 
when fy (~)=mz, | 


пе 
( Z4 dT м2 (2,7°-+2207+... 222, | £a) 
(n—1) 
= | n ue (А MET r) (А14) 
0 2 a 5n 
and 
ne (пе (n—1)e  (n—1)e 
i f 2207 дт |, | 2207 07 
0 Jo 0 0 
€ (n—1)e пе | 
s( f ЕВ |, гаг) (А15) 


An improved value for the velocity is obtained from the 
expression 


теа, Ро (2а, a Es) (A16) 
This value is used in the оов of the estimated velocity 
and displacement 25 ,, for the next interval. - 

If desired, improved values of the displacement and 
acceleration for the nth interval subsequent to the beginning 
of shock-strut deflection can be obtained as follows: 


€ А А 
Zay T 22.15 Со р 22) 


2 
га, Ней, а (22,_,+22,) (A17) 


and 


= (22, 22, та) 


where (22, 22,, Та) is ап appropriate equation for the system, 
such as equation (A7). 

With the values of 2, 2;,, and 2 „ the motion variables for 
the upper mass Z,,, 2, and 2, can be calculated separately 
from equations (A4), (A5), and (A6), as previously noted. 

То setting up the numerical procedure used in obtaining the 
solutions presented in this report, an evaluation of the errors 
introduced by the procedure indicated that it would not be 
necessary to calculate the improved values of the displace- 
ment 2, (eq. (A17)) or the acceleration Z, (eq. (A18)). How- 
ever, improved values of the velocity %, were calculated by 


(A18) 
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means of equation (А16) for the purpose of determining 
estimated values of the velocity 2 and the displacement ж 
(eqs. (А11) and (А12)) for the increment immediately 
following. 

In order to illustrate the application of the method, a 
tabular computing procedure for the solution of the system 
represented by equations (А1), (A2), and (АЗ) is presented in 
table I. 

QUADRATIC PROCEDURE 

In this step-by-step procedure a quadratic variation of 
displacement is assumed over successive equal finite time 
intervals for the purpose of extrapolating values of the 
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motion variables from one interval to the next. With this 
assumption the displacement variation over two successive 
equal time intervals is completely determined by the three 
values of displacement at the beginning and end of each of 
the two intervals. By writing the quadratic variation in 
difference form, the velocity and acceleration at the midpoint 
of the double interval can be expressed in terms of the three 
displacement values previously mentioned. Substituting for 
the velocity and acceleration in the differential equations for 
the system yields difference equations of motion in terms of 
successive displacement values which can be evaluated 
interval by interval. 


TABLE I 
LINEAR PROCEDURE 


Row Quantity 


e Zin 
© Ву (аа) | | __ 202 
(o |, Рт, (4+ 
© | p Pra dr 

| 해 
9 z 
®© 22, 

1 Е 
2 (23, t зал) 
(9) Zu 
AL 

© z 
@ 2. 
3 Zi. 


ТО, denotes value for previous time interval. 


Equation 


Equation (A14) 
Equation (A15) 
Equation (A7) - 
fags tS (аа EON 
b, e 8,) 
за (ena th) 


Equation (A4) 


Equation (A5) 


Equation (А6) 


Procedure f 


(2,- €; 


Determined from tire force- 
deflection characteristics. 


954-5 [0+0 
Given by equation (A7). 
О, +: [0+0 
@+51©+@sle 


Q--5 Ф+ Фе 


Given by equation (A4). 


Given by equation (A5). 


_ Given by equation (Аб). 
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The following derivation shows how the procedure can be 
applied to the determination of the behavior of the landing 
gear subsequent to the beginning of shock-strut deflection at 
time Е. In order to simplify the derivation, internal friction 
forces are again neglected in setting-up the equations of 
motion. 

The assumption of a quadratic variation of displacement 
with time (constant acceleration) over two successive inter- 
vals, each of duration e, permits expressing the velocity and 
acceleration at the midpoints of the double interval (see 
sketch) in terms of the displacement values at the beginning, 
midpoint, and end of the double interval by the equations 
(see ref. 5, p. 16): 


€ 


a таз (A19) 
and 
Zuni — 22, 2. af 
ccs я--1 лаж n—1 (A20) TuÍ-f, 


where 2,, 2,, and 2, are the velocity, acceleration, and displacement at the end of the nth interval (т--те) after the 
beginning of shock-strut deformation and z,., and 2,4, are the displacements at the end of intervals n—1 and n-+1, 
respectively. 

Substituting the difference relations for 21, 2, 21, and % into equations (А1) and (АЗ) permits writing the equations 
of motion for the landing gear in difference form as follows: 


ze (21,4124, 71,..) аа £1, 4 22, ға, 1) B [17 0(21,—23,)] "+ D—0 (A21) 
and | 
21444177 221 21, 4— 0 (22,4, — 222, + 27000 -Не| Ру.(2з) +E) . (A22) 


where the constants are as defined in the previous section. 


Substituting for 21 ,, in equation (А21) gives 


я 1 | 
ааа. + РА WwW? | wem. —44Wy?W«(g.A ИВ. t WW) —gA W (AW oni ВЕ Ау | | (А28) 


where 
Опа = 2 W323, — 22.1 8e | Е (2з,) TE] 


Ва-ызт 20725, +(Wi— Wa)Z2,_, + 2Wi(zi,— 21,21) —ge"[ Py, (22,) TE] 


and 


2.2 
ㆍ veu ETE. [Bp O(21, —23,)] FD} 


Equations (А22) and (A23) are essentially extrapolation formulas which permit the determination of values for the 
upper-mass and lower-mass displacements to come from the values of displacement already calculated. These equations 
thus permit step-by-step calculation of the displacements as the impact progresses, starting with the initial conditions, 
from which the upper-mass and lower-mass velocities and accelerations can be determined by means of equations (A19) 
and (А20). | 

Since the calculation of the displacements г and z at any instant by means of equations (422) and (А23) requires 
values for the displacements at two previous instants, the routine application of these equations can begin only at the 
end of the second interval (r—2e) following the beginning of shock-strut deflection. Before the displacements at the end 
of the second interval can be calculated, however, it is necessary to determine the displacements at the end of the first 
interval. These values can be obtained from the conditions of motion at the instant of initial shock-strut deflection by 
applying equations (А19) and (А20) to the instant і. 
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At the instant of initial shodk-strut deflection . 


"loud 295.0 Zr 


(A24) 


Application of the difference equations (A19) and (A20) 
to the instant і--і, (that is, n=0) gives the following 
equations: 


Znwl' "Сбла--1 
дє 


та 22--Р Z aw—l 


5 Znel 
s 
«2 


i= 
(A25) 


Since the landing gear is considered as а one-degree-of- 
freedom system from initial contact up to the instant і--і;, 
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the foregoing application of the difference equations results 
in identical values for the upper-mass displacement and 
lower-mass displacement at the end of the first interval. 
Simultaneous solution of equations (A25) gives the following 


expression for the displacement at the end of the first interval: 


(А26) 


ааа 5 


With the values for z and z,.:, equations (А22) and 


. (А28) permit the step-by-step calculation of the upper-mass 


and lower-mass displacements subsequent to the first interval 
following the beginning of shock-strut deflection. The 
corresponding velocities and accelerations of the upper and 
lower masses can be determined from the calculated displace- 
ments by means of equations (A19) and (A20), as previously 
noted. 

A tabular computing procedure illustrating the application 
of the method is presented in table II. 


TABLE II 
QUADRATIC PROCEDURE 


Quantity 


Тоне ed 개다 재개 매매 192 ed 


am ат что еи "з= чыл AXE UPS ӘР UNA шы Wm ч== ==” чыз 


— «шу «кы. TI ee ll бын қыз ана а | 


a 


Equation (А23) 


Equation (А22) 


— 2z; Жан 21 


Procedure { 


LIP 으니 ee ee ee ee 고 으 에 호 ee 호 은 ФМИ «ер; өше че Деби LENIN E 


— G, 


… | Ф, 


Given by equation (A23). 


Given by equation (А22). 


м--1 


1 © denotes value for previous time interval. 
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RUNGE-KUTTA PROCEDURE 


In this step-by-step procedure the differences in the de- 
pendent variables over any given interval of the independent 
variable are calculated from a definite set of formulas, the 
same set of formulas being used for all increments. ‘Thus the 
values of the variables at the end оѓ any given interval ате 
completely determined by the values at the 600 of the рге- 
ceding interval. Unfortunately, however, unless the equa- 
tions to be integrated are relatively simple, the method can 
become quite lengthy. 

The following derivation illustrates the application of the 
Runge-Kutta method to the generalized equations of motion 


(eqs. (21)) for the simplified system considered in the section’ 


on generalized results. Since these equations can be readily 
reduced to the first order, they can be integrated by the 
step-by-step application of the general equations given on 
pages 301 and 302 of reference 6 for first-order simultaneous 
differential equations. 

The generalized equations for the simplified system pre- 
viously discussed (eqs. (21)) are 


(uy! — и)" =0 
(u;’ —u y — w= 0 


wu"! Hua = 0 


Inasmuch as any two of these equations are sufficient to 
describe the behavior of the system, only the last two equa- 
tions are employed in this procedure. These equations can 
be reduced to a first-order system by introducing the new 


variable i 
W= U (А27) 
so that 
W =u," (A28) 
and the equations of motion become 
(w» —us)—u,—0 
| (А29) 
w’-+m=0 


Solving equations (А29) for из” апа ш”, respectively, gives 
и! = 0 — ид (A30) 
Vp! = — U | (A31) 


Applying the general procedure presented in the reference 
previously cited to the simultaneous equations (À27), (À30), 
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and (À31) gives 
Аш, mm (kit Bhat hy ES 
Аара w, ае (ht let Qa tly) (A32) 


| ‚1 
Ata —U2, 176 (m, 4- 2mm34-2m3-4- m.) 
where 
bh —7w. A0 
| l 
la (пила) АВ 
2 
CONSE) А0 


ky=(Wy_-1+ 60 А0 


(тш — Uz, „Ад 
m 
ЕСТЕ A0 


ш--(ш, а) A0 


l= — (ua, .,- ma) A0 


m= (Юта ү) Абиз" А6 
ma=| (ve) аи, + ав 
СТО NEC | 40 


т. | (0. 14-45) — из, ms | Ад 


With this procedure, ш, w, aud uz can be calculated in 
step-by-step fashion from the values for the preceding inter- 
val, the procedure beginning with the initial conditions. 
From these values, w’, «^, and и” can be calculated by 
means of equations (A27), (A28), and (A30), respectively. 


--- m- — -— ——— 
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APPENDIX B 


SOURCE OF EXPERIMENTAL DATA 


Following is a brief description of the apparatus and test 
specimen used in obtaining the experimental data presented 
in this report. 

EQUIPMENT 

The basic piece of equipment employed in the tests is the 
carriage of the Langley impact basin (ref. 7) which provides 
means for effecting the controlled descent of the test speci- 
men. In these tests.the impact-basin carriage was used in 
much the same manner аз а, conventional stationary landing- 
gear test jig (see ref. 8). In order to simulate mechanically 
the wing lift forces which sustain an airplane during landing 
the pneumatic cylinder and cam system incorporated in the 
carriage was used to apply a constant lift force to the dropping 
mass and landing gear during impact. Тһе lift force in these 
tests was equal to the total dropping weight of 2,542 pounds. 


TEST SPECIMEN 


The landing gear used in the tests was originally designed 


for a small military training airplane having a gross weight of 
approximately 5,000 pounds. The gear is of conventional 
cantilever construction and incorporates a standard type of 
oleo-pneumatic shock strut. The wheel is fitted with a 27- 
inch type I (smooth-contour) tire, inflated to 32 pounds per 
square inch. The weight of the landing gear is 150 pounds. 
The weight of the lower mass (unsprung weight) is 131 
pounds. 

In the present investigation the gear was somewhat modi- 
fied in that the metering pin was removed and the original 
orifice plate was replaced with one having a smaller orifice 
diameter. Figure 14 shows the internal arrangement of the 
shock strut and presents details of the orifice. Other perti- 
nent dimensions are presented in table ПП. The strut was 
filled with specification АМ-УУ-0-366В hydraulic fluid. 
The inflation pressure with the strut fully extended was 43.5 
pounds per square inch. In these tests the landing gear was 
mounted with the shock-strut axis vertical. Figure 15 
is a photograph of the landing gear installed for testing. 


TABLE ПІ 


IMPORTANT CHARACTERISTICS OF LANDING GEAR 
USED IN TESTS 


AEE 기라 КО наннан OES 0. 04708 
「 Tr Same ER AOI A RE Ne CES 0. 0005585 
А Tere Seto ee ae ee eat 0. 03545 
ТОРН E RE сын сыз Ne БЕР 6, 264 
Пе EE OEE a КЕТЕТИН A SE eR Ie 0. 5521 
Газ (eee HEP OL dot E UNS 2. 22604 
И Ic c oin iltud il 2,411 
гид | «a мс Ue АР ы QUIS DAMM DAN 131 


INSTRUMENTATION 


A variety of time-history instrumentation was used during 
the tests. The vertical acceleration of the upper mass was 
measured by means of an oil-damped electrical strain-gage 
accelerometer having a range of --8g and a natural frequency 
of 85 cycles per second. A low-frequency (16.5 cycles per 
second) NACA air-damped optical-recording accelerometer, 
having a range of —ig to 6g, was used ав а stand-by instru- 
ment and as a check against the strain-gage accelerometer. 
Another oil-damped strain-gage accelerometer, having a 
range of +12¢ and а natural frequency of 260 cycles per 
second, was used to determine the vertical acceleration of the 
lower mass. The vertical displacement of the lower mass 
(tire deflection) and the shock-strut stroke were measured 
separately by means of variable-resistance slide-wire poten- 
tiometers. "The vertical displacement of the upper mass was 
determined by addition of the strut-stroke and tire-deflection 
measurements. The vertical velocity of the landing gear at 


theinstant ofground contact was determined from the output 


of an elemental electromagnetic voltage generator. A time 
history of the vertical velocity of the upper mass was ob- 


‘tained by mechanically integrating the vertical acceleration 


of the upper mass subsequent to the instant of ground con- 
tact. Electrical differentiation of the current output of the 
strut-stroke circuit provided time-history measurements 
of the shock-strut telescoping velocity. The instant of 
ground contact was determined by means of а micro- 
switch, recessed into the ground platform, which closed 
a circuit as long as the tire was in contact with the platform. 

The electrical output of the instruments was recorded on a 
14-channel oscillograph. The galvanometers were damped 
to approximately 0.7 critical damping and had natural fre- 
quencies high enough to produce virtually uniform response 
up to frequencies commensurate with those of the measuring 
instrumentation. A typical oscillograph record is shown in 
figure 16. 

It is believed that the measurements obtained in the tests 
are accurate within the following limits: 


Measurement Accuracy 
Upper-mass acceleration, 2-.-----.-.--.---.---..-........ +0. 2 
Force on upper mass, Б-------------<---<-<-<-<-<-< +600 
Lower-mass acceleration, Е ----------------- P +0. 8 
Vertical velocity at ground contact, Їрз------------------- +0. 1 
Upper-mass velocity during impact, Їр8-.-.--------------- +0. 5 
Upper-mass displacement, В .-------------------------- +0. 05 
Lower-mass displacement, Ї{_--.—.----------------------- +0. 03 
Shock-strut stroke, f£... LLL LLL LL LL c LLL Lecce ceu +0. 03 
Shock-strut telescoping velocity, Ёр3-.-------------------- +0. 6 
Time alter Contact, 866.2 ou co cea eis eid nm dae mas eia ae es +0. 003 
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Етсове 14.—Shock strut of landing gear tested at Langley impact basin. 
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FrcunE 15.—View of landing gear and instrumentation. 
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Ficure 16.—Typical oscillograph record obtained during test in Langley impact basin. 
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